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In this paper, we compare the relative sensitivities of gamma-ray and neutrino observations to the dark mat- 
ter annihilation cross section in leptophilic models such as have been designed to explain PAMELA data. We 
investigate whether the high energy neutrino telescope IceCube will be competitive with current and upcom- 
ing searches by gamma-ray telescopes, such as the Atmospheric Cerenkov Telescopes (HESS, VERITAS and 
MAGIC), or the Fermi Gamma-Ray Space Telescope, in detecting or constraining dark matter particles anni- 
hilating in dwarf spheroidal galaxies. We find that after 10 years of observation of the most promising nearby 
dwarfs, IceCube will have sensitivity comparable to the current sensitivity of gamma-ray telescopes only for 
very heavy (mx ^ 7 TeV) or relatively light (mx ^ 200 GeV) dark matter particles which annihilate primar- 
ily to p + fi~ . If dark matter particles annihilate primarily to t + t~, IceCube will have superior sensitivity only 
for dark matter particle masses below the 200 GeV threshold of current Atmospheric Cerenkov Telescopes. 
If dark matter annihilations proceed directly to neutrino-antineutrino pairs a substantial fraction of the time, 
IceCube will be competitive with gamma-ray telescopes for a much wider range of dark matter masses. 

PACS numbers: 95.35.+d;95.30.Cq,98.52.Wz,95.55.Ka UTTG-13-09 TCC-030-09 FERMILAB -PUB -09-5 89-A 



INTRODUCTION 

Of the many dark matter candidates to have been proposed, those within the classification of Weakly Interacting Massive 
Particles (WIMPs) are often considered to be the best motivated (for reviews, see Refs. 0] El)- In particular, if a stable par- 
ticle with an electroweak-scale mass and electroweak-scale interactions exists, such particles would have annihilated among 
themselves at a rate in the early universe that naturally led them to freeze-out with a relic density similar to the measured dark 
matter abundance. This process of dark matter annihilation is also predicted to be taking place in the present universe, providing 
the basis for dark matter indirect detection experiments. Such experiments search for the annihilation products of dark matter 
particles, including electrons and/or positrons, antiprotons, photons, and neutrinos. Promising sites for the observation of dark 
matter annihilation products include the core of the Sun Q, the Earth [4 |, our Galactic halo [5 |, Galactic center [6 |, and dwarf 
satellite galaxies (7). Here, we examine the sensitivity of the IceCube neutrino detector to dark matter annihilations taking place 
in dwarf satellite galaxies in the Milky Way. 

In the case of dark matter candidates that annihilate primarily to gauge bosons or hadronic final states, gamma-ray telescopes 
provide a more sensitive test of dark matter annihilations taking place in dwarf spheroidals than can be accomplished with 
existing or planned neutrino telescopes. This is not necessarily the case, however, if the dark matter annihilates largely to 
leptons. Dark matter annihilating to leptonic final states has become increasingly well motivated by the anomalous observations 
recently reported by several cosmic ray experiments, including PAMELA, which has observed an excess of cosmic ray positrons 
(relative to electrons) between 10 and 100 GeV lfl4l (along with previous experiments, including HEAT 1 15| and AMS-01 1 16|, 
which also reported evidence for such an excess). A surplus of cosmic ray electrons and/or positrons has also been reported by 
the Fermi Gamma-Ray Space Telescope (FGST) [17], however this excess is less than that previously reported by ATIC lfl8ll . 

Although an excess of cosmic ray positrons/electrons may come from more conventional sources such as pulsars lfT9l , a great 
deal of interest has been generated in the possibility that these signals might result from dark matter particles annihilating in 
the local halo of the Milky Way. Efforts to produce such signals with dark matter, however, have faced some model-building 
challenges. First, as there is no observed excess of cosmic ray antiprotons, the annihilation channels must be "leptophilic," 
i.e. the WIMPs must annihilate preferentially to leptons. Many models have been proposed with this property |20l . Second, 
the local halo density of dark matter is insufficient to explain the signals unless the annihilation rate is supplemented by a large 
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factor ~ 10 — 10 3 relative to that predicted for a typical smoothly distributed thermal relic. Such a boost may arise due to the 
properties of the dark matter candidate itself; for example, via a nonperturbative Sommerfeld enhancement to the low-velocity 
cross section resulting from the exchange of a light state ETI . or a Breit-Wigner enhancement [22, 23 1. Because of the smaller 
velocity dispersion of dark matter particles in smaller halos, the annihilation cross section in a dwarf galaxy may be even higher 
than that in the halo of the Milky Way fl2Tl 12414261 . Alternatively, an enhancement to the annihilation rate over the standard 
prediction can result if the dark matter is not distributed smoothly. In this paper, we examine bounds on the observed WIMP 
annihilation cross section, (crv), including any relevant enhancements or boost factors, the origin of which we leave unspecified. 

The prospects for IceCube observations of high energy neutrinos from dark matter annihilations or decays in the Galactic 
halo have been considered in [8-11]. It was found that IceCube will be sensitive to WIMP annihilation cross sections of order 
10 -24 to 10 -23 cm 3 s _1 , depending on the WIMP mass. As Super-Kamiokande is located in the Northern Hemisphere, and 
therefore is subject to significantly reduced atmospheric backgrounds from the direction of the Galactic center, its sensitivity to 
the neutrino flux from dark matter annihilations in the inner Milky Way has been studied in |[T2"1 [T3l and found to be roughly 
an order of magnitude less sensitive than the IceCube projections. Nonetheless, Ref. iPPUl finds that Super-Kamiokande already 
constrains annihilations to t + t~ as the explanation for the cosmic ray anomalies (though this channel is also disfavored by 
the large expected 7 flux), while Ref. |fl3l finds that Super-Kamiokande may also constrain the scenario in which neutrinos are 
directly produced in the annihilations as frequently as charged leptons. While these studies are concerned with the signature of 
annihilations occurring in the Galactic center, here we focus instead on the annihilations taking place in satellite dwarf spheroidal 
galaxies of the Milky Way. In particular, we estimate the neutrino fluxes from WIMP annihilation in the dwarf galaxies Draco, 
Willman 1, and Segue 1, and examine the prospects for their observation with the IceCube neutrino telescope located at the 
South Pole. 



DARK MATTER ANNIHILATION IN DWARF SPHEROIDAL GALAXIES 



In general, the differential flux of particles of type j from the annihilation of dark matter particles of mass mx is given by 

dEj 8nm x ^ Jt dE j>F y ' ' 

where fp is the fraction of annihilations which produce a final state F, and dNjp/ dEj p is the differential spectrum of particles 
j from an annihilation to final state F. J(Afl) is the square of the dark matter density integrated along the line-of-sight, averaged 
over the solid angle AO: 



J(Afi) = / dSl f Px (s)ds. (2) 

A " ./AO Jl.o.s. 



Dwarf spheroidal galaxies are promising sources for indirect dark matter searches for a number of reasons. First, they con- 
tain relatively large dark matter densities, and thus may produce sizable fluxes of annihilation products. Furthermore, dwarf 
spheroidals contain relatively little in the way of baryonic material (stars, gas, etc.) and are astrophysically simple. The typical 
mass of known dwarf spheroidal galaxies is ~ 10 7 M , which is distributed over a volume on the order of a cubic kiloparsec. 
Dwarfs are largely devoid of astrophysical activity, and have very large mass to light ratios, Segue 1 being a particularly extreme 
example l27l . The systems therefore lack astrophysical sources which could potentially mimic a signal of dark matter annihila- 
tion. Although simulations suggest that many more are likely to exist, approximately 25 dwarf galaxies within the local group 
have been discovered thus far. 

The dark matter distribution p(r) within a dwarf spheroid can be fit with a five parameter density profile ll28l 

P(r) = ^ ^r, (3) 

(Fl) a (l + (^) b )^ 

where r is the distance from the dynamical center of the dwarf galaxy, r s is the scale radius, and po is the central core density. 
Typical ranges for the parameters a, c, and b, which determine the inner slope, outer slope, and transition between the two, 
respectively, are: a = [0.0 — 1.5], b = [0.5 — 1.5], and c = [2 — 5]. A^-body simulations find cusped inner profiles with 
a = [1.0 — 1.5]. The profile of each individual halo can vary, however, depending upon its own merger history. Given this 
uncertainty, we follow Ref. l28l which marginalized over these five parameters, including r s and po, over the ranges specified 
above with flat priors. Although the parameter ranges are generous, Ref. ESI finds that the velocity dispersion data fix the 
line-of-sight integral, J(Af2), to lie in a relatively small range. For details, see Ref. [29 1. 

Presently, more than half a dozen dwarf galaxies are known which could potentially provide an observable flux of neutrinos. 
These dwarfs are relatively near the Solar System (tens of kpc) and, importantly, are in the northern hemisphere. This enables 
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IceCube's background of atmospheric muons to be largely avoided by looking through the Earth, towards the northern hemi- 
sphere. We have selected three nearby dwarfs as sources of interest for IceCube; Draco, Willman 1, and Segue 1. In Table[I] we 
display the relevant properties of these satellite galaxies l30l . with J(Af2) Ail and errors representing the 2a range taken from 
Ref. ||3T| . We note that there is some uncertainty in the observational status of Willman 1; it is currently unclear whether it is a 

TABLE I: Properties of selected Milky Way dwarf galaxies. 
Galaxy Distance (kpc) Logio J(Afi) Afl Declination 
Draco 80 18.63 ± 0.73 +57°54'55" 

Willman 1 38 19.55 ±1.19 +51°03'00" 

Segue 1 23 19.88 ± 0.82 +16°04'55" 

dwarf galaxy or a globular cluster. It follows that there is considerably more uncertainty in the value of J (Ail) Ail. 

IceCube is designed to have an angular resolution of approximately 1° for muon tracks 1321 . Given this resolution, the dwarf 
galaxies in Tableware sufficiently small and distant that we can treat them as point sources. Atmospheric Cerenkov Telescopes 
and the Fermi Gamma Ray Space Telescope have angular resolutions much smaller than IceCube. For example, VERITAS has 
observed Draco, Ursa Minor, and Willman 1 over regions within 0.15° of each galaxy's center 11331 . 

Several dwarf galaxies have been observed by ACTs as potential sources of gamma-rays: Draco by MAGIC [34| and 
STACEE [35 1, Sagittarius (which is visible from the Southern Hemisphere) by H.E.S.S. [36 1, and Draco, Willman 1, and Ursa 
Minor by VERITAS [ 33 ]. None of these telescopes have observed any significant signal and, therefore, have placed upper limits 
on the flux of gamma-rays coming from these sources. The strongest constraints on the annihilation rate of WIMPs in dwarf 
galaxies come from H.E.S.S. and VERITAS [28]. H.E.S.S. observations of the center of the Sagittarius dwarf yield an upper 
limit on the gamma-ray flux above 250 GeV of $(i? 7 > 250 GeV) < 3.6 x 10~ 12 cm _2 s _1 at the 95% confidence level 031 . 
VERITAS observations of Draco, Ursa Minor, and Willman 1 bound the gamma-ray flux from each of these objects above 
200 GeV to be $(£ 7 > 200 GeV) < 2.4 x 10~ 12 cm^s" 1 11331. 



MUON NEUTRINOS FROM DWARF SPHEROIDAL GALAXIES 

In this study, we consider three leptophilic scenarios: WIMPs which annihilate solely to t + t~, solely or to both 

and v^v^, each with branching fractions of 50%. The differential flux of neutrinos from annihilations directly to neutrino- 
antineutrino pairs in a dwarf galaxy is given by Eq. [TJ and the differential spectrum at production takes the very simple 
form 

dN v dN^ 

— — = — —=6{E V -m x ). (4) 
dE Ufi dE^ V X> 

In the cases in which the WIMPs annihilate to t + t~ or /i + /J~, we use PYTHIA ||37l to calculate the resulting neutrino spectra. 
In all cases, we include the effects of three-flavor vacuum oscillations [38 1. 

Atmospheric muon neutrinos constitute the most serious background for detection of from dark matter annihilation prod- 
ucts from northern hemisphere dwarf spheroidal galaxies. This background is a function of the zenith angle observed |39|. We 
use Ref. l40l to calculate the atmospheric neutrino background specific to the direction of each of the three source dwarfs. We 
note that the predicted background spectra are in agreement with the observations of Amanda-II [41 1. 

The rate of muon tracks from charged current neutrino interactions observed at IceCube is calculated by combining the in- 
coming spectrum of muon neutrinos with the probability of those neutrinos being converted to muons above the energy threshold 
of the telescope [42], which is given by 

P(E V ,E^% = p N A / dy^(E„, y) [R^E V (\ - y), E^) + L]. (5) 
Jo a V 

Here, p = p- lcc ss 0.9 g cm~ 3 is the density of the target medium (ice) and Na = 6.022 x 10 23 is Avogadro's number. The total 
charged current cross section for neutrino nucleon scattering, ovjv, can be taken as 

Vvn(E v ) + (T„n{E v ) _ ne ^ in -36 ( E u \ 096 2 



3 - 06x10 loitvj cm (6) 

for 100 GeV < E v < 1 TeV, and is approximately flatly distributed in y |43l . Finally, + L is the effective detector size, the 
sum of the physical length of the detector, L, and the distance a muon travels before its energy falls below the threshold of the 
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experiment, R^. The muon range in ice is 

2.0 + 4.2(^/TcV) 



R,j, m 2.4 km x In 



_2.0 + 4.2(££ hr /TeV) 
The event rate for both the background and the signal is given by the integral 

N = J A cS P{E v ,E^)^dE v . 



(7) 



(8) 



The effective area, A e s , in IceCube is approximately one square kilometer. The DeepCore supplement to the IceCube detector is 
designed to increase the sensitivity to low energy neutrino-induced muons with 10 GeV < E^ < 100 GeV, with the improvement 
most substantial at the lowest energies B4ll . However, the dark matter annihilation signal in IceCube and DeepCore is dominated 
by muons with energies close to the dark matter mass: for the lowest WIMP masses considered here, mx — 100 GeV, the signal 
comes primarily from neutrinos with 50 GeV < E v < 100 GeV. In this case, considering only the IceCube and DeepCore 
sensitivities reported in Fig. 4 of Ref. l44l . we expect an improvement of roughly 20% over IceCube alone. 

WIMP masses below ~ 100 GeV are not compatible with the anomalous cosmic ray excesses. However, the sensitivity to 
neutrinos with E v ~ 10 GeV is improved by including the DeepCore addition by more than an order of magnitude over that of 
IceCube alone (44). At the same time, since the maximal angle between a muon track and the primary neutrino increases at low 
neutrino energies as l/y/E v , the angular resolution is significantly degraded at low energies, resulting in a decreased signal to 
background ratio for localized objects such as the dwarf satellite galaxies in Table [I] We therefore estimate that the DeepCore 
addition would improve the sensitivity to low-mass dark matter annihilating in dwarf satellite galaxies by, at most, a factor of 
a few. However, as the cosmic ray anomalies do not favor light leptophilic dark matter, we report only on the sensitivity of the 
IceCube neutrino detector to dark matter annihilations in dwarf galaxies for WIMPs with mx > 100 GeV. 

In comparing the event rate to that predicted from the background of atmospheric neutrinos, we include events with energy 
up to 25% greater than the WIMP mass, which conservatively accounts for the finite energy resolution of IceCube (32). We 
can now calculate the exposure time, T, necessary to achieve a detection with a given level of statistical significance over the 
background of atmospheric neutrinos. Since the background is essentially isotropic over the angular window of interest we 
follow the analysis of Ref. fl45l to find the minimum exposure time, 



■T 1 ^atm + ^DMA „ 

r = <w — ^ > <- 9 ) 

DMA 



where Ndma is the event rate from dark matter annihilations in the dwarf and N atm is the event rate from the atmospheric 
background. In the following section, we use this approach to determine the sensitivity of IceCube to dark matter annihilations 
taking place in the dwarf galaxies Draco, Willman 1, and Segue 1. 



PROSPECTS FOR ICECUBE AND COMPARISON TO GAMMA-RAY SEARCHES 

We now compare the relative merits of gamma-ray and neutrino observations in their ability to place limits on the dark 
matter annihilation cross section in leptophilic models. In Fig. [T] we present IceCube's sensitivity to neutrinos from dark 
matter annihilations in Draco, after 10 years of observation. In the left panel, we show results for the annihilation channels 
XX — > /i + /J~ in red (lower contours) and XX — > t + t~ in blue (upper contours). As a consequence of neutrino oscillations, 
red and blue contours are distinct only at low mx- The thick dashed lines denote the projected 2a upper limit on the annihilation 
cross section as a function of mass (using the 2a low value for J(AO)AQ), whereas the shaded regions bounded by thin solid 
lines denote the minimal cross sections necessary for a 5cr discovery, given the ±2c range of values for J(AO) Af2. In the right 
panel, results are shown for the case in which XX —> fi + n~ and XX — > v^v^ each occur with branching fractions of 50%. 
One can see that, due to the uncertainty in the dark matter halo properties, the potential discovery regions extend well below the 
exclusion contours. 

For comparison, in each frame we also show the corresponding upper limits on the dark matter annihilation cross section 
from VERITAS observations of Draco, as found in Ref. (28). In the left panel, the grey and black dotted curves represent the 
upper limits on (av) assuming annihilations to /J + /i and t + t~, respectively, while in the right panel, the grey dotted curve is 
the VERITAS upper limit assuming XX — > /J + /i~ and XX — > v^v^ with equal branching fractions. We note that Ref. ||28) 
uses the 90% confidence level lower limit on J(Af2)Af2, while we consider the slightly more generous 2a lower limit. Our 
projected sensitivities are therefore slightly more conservative. As Draco is the least bright of the dwarf galaxies considered 
here, we find that for WIMP masses > 200 GeV, VERITAS has already constrained the annihilation cross section roughly as 
tightly as IceCube will with 10 years of data for annihilation purely to /i + fi~ or t + t~. In the case of annihilations to 
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FIG. 1: Projected upper limits on the dark matter annihilation cross section, (cv) (alternatively, boost factor), at 95% confidence level as a 
function of mass (thick dashed) from neutrino observations of the dwarf spheroidal galaxy Draco. Also shown are the 5o discovery regions 
(shaded). In the left frame, the red (lower) region and lines denote the case in which dark matter annihilates to p + p~ , whereas the blue 
(upper) regions and lines denote annihilations to t + t~ . The right frame is for the case in which dark matter annihilates to both fi + n~ and 
v^Vfi with equal probability. In each frame, the dotted lines denote the current limits from gamma-ray observations by VERITAS. See text for 
more details. 
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FIG. 2: The same as Fig. 1, but for observations of the dwarf spheroidal galaxy Willman 1. 



IceCube may place competitive limits on the WIMP annihilation cross section only for mx Si 5 TeV or mx 200 GeV. For 
annihilation to t + t~ the current constraints on (av) from VERITAS observations of Draco are stronger than any IceCube will 
be able to achieve in the next 10 years for mx > 200 GeV. This stems from the fact that taus often decay hadronically, resulting 
in many more high energy photons than in the purely leptonic decays of muons. For annihilation to taus, a substantial portion of 
IceCube's discovery reach has already been excluded by VERITAS. 

Instead of presenting our results as limits (or discovery reach) on the dark matter annihilation cross section, we can alterna- 
tively describe them in terms of the boost factor to the annihilation rate in a given dwarf spheroidal galaxy. Along the right side 
of each frame, the constraints and discovery prospects are given in terms of this quantity, which we define as the annihilation rate 
divided by that predicted for dark matter with {av) = 3x 10~ 26 cm 3 /s and distributed smoothly over the halo. This definition 
allows us to leave the origin of the boost factor unspecified. 

Turning to the right panel, we find that if neutrinos and muons are each are directly produced in 50% of annihilations, IceCube 
will have comparable or better sensitivity than VERITAS for all WIMP masses. In fact, in this case IceCube's discovery reach 
is largely untested by ACTs. We note that if the branching fraction to neutrinos is larger than 50%, IceCube's sensitivity will be 
further increased, while that of ACTs and other gamma-ray telescopes will be reduced. 

In Fig. [2] we present the same information as Fig.[T] but for the dwarf galaxy Willman 1. Willman 1 is somewhat brighter than 
Draco, and therefore we expect all sensitivities to improve. However, there is also considerably more uncertainty in the halo 
profile, leading to large discovery regions and inflated prospective limits. For annihilations to charged leptons only, boost factors 
of O(10 2 ) may be accessible at IceCube if J(Afi)Af2 is near its 2a upper limit. For low mx, these annihilation cross sections 
have not yet been probed by ACTs. For annihilation to and v^v^, shown in the right panel, boost factors as low as O(10) 

may be accessible in the most optimistic case. 
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FIG. 3: The same as Fig. 1, but for observations of the dwarf spheroidal galaxy Segue 1. In this figure, the dotted lines denote the preliminary 
limits from the Fermi Gamma-Ray Space Telescope A46V . 



Finally, we examine the closest of the three dwarfs we have considered, Segue 1. In Fig. [3] we present IceCube's projected 
sensitivity to dark matter annihilations in Segue 1, along with the current upper limits on the annihilation cross section from 
FGST measurements l46l . For annihilation to muons (taus), IceCube may eventually set stronger exclusion limits for mx > 
450(750) GeV. If there is a 50% branching fraction to neutrinos, IceCube will do better for all WIMP masses shown. As Segue 
1 is newly discovered, we expect that ACTs will observe this galaxy in the near future. 



DISCUSSION AND CONCLUSIONS 



In this article, we have considered the prospects for the kilometer-scale neutrino telescope IceCube to detect neutrinos from 
dark matter annihilating in dwarf spheroidal galaxies, and have compared these to the constraints placed by gamma-ray tele- 
scopes. We find that if dark matter annihilates primarily to muons, taus, and/or neutrinos, IceCube can potentially provide 
constraints comparable to or even stronger than those obtained by current gamma-ray telescopes. 

Other dwarf spheroidal galaxies, such as Ursa Minor, which yield weaker gamma-ray limits on the dark matter annihilation 
cross section will also be less promising to be observed by IceCube. An analysis of the neutrino flux from all dwarf galaxies 
accessible to IceCube could potentially improve the sensitivity to the annihilation cross section, however, the maximal improve- 
ment assuming all dwarfs are equal is proportional to 1/ y/N, where N is the number of dwarfs in the sample. All dwarfs, of 
course, are not equal. In order to carry out such an analysis, one must assume that the mechanism responsible for the boosted 
(relative to thermal) annihilation cross section results in the same observed cross section in each dwarf galaxy. As the velocity 
dispersion in each dwarf is independent and the dumpiness of each dwarf is unknown, there is no reason to expect that the 
observed annihilation cross section in any two dwarfs should be the same, thus, at this point, the only appropriately model- 
independent interpretation is to view the annihilation signal from each dwarf independently. Given the velocity dispersion in 
each dwarf, one could derive limits on particular dark matter models assuming a specific velocity-dependent enhancement to the 
annihilation cross section. Such a study may be useful in the near future. 
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